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O ' Abstract 

' This paper deals with existence and uniqueness, in viscosity sense, of a solution for a system of m 

variational partial differential inequalities with inter-connected obstacles. A particular case of this 
system is the deterministic version of the Verification Theorem of the Markovian optimal m-states 
switching problem. The switching cost functions are arbitrary. This problem is connected with the 
valuation of a power plant in the energy market. The main tool is the notion of systems of reflected 
backward stochastic differential equations with oblique reflection. 
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1 Introduction 

O 

The multi-modes switching problem is by now well documented both in the economics or mathematical 
literatures (see e.g. [H El [Q El HOI US [HI [El [2Ql [271 ISSl [29], etc. and the references therein). The 
pioneering work of Brennan and Schwarz [4J deals with a two-modes switching problem describing the 
life cycle of an investment in the natural resource industry. A major switching problem of interest 
is related to the energy market. Actually let us consider a power plant which has several modes of 
production and which is put in a specific mode according to its profitability which depends on the 
electricity price in the market. The manager of the plant aims at maximizing its global profit. For this 
objective, she implements an optimal strategy S* which is a pair of two sequences (rk)k>i and (£fc)jfc>i 
describing respectively the optimal successive switching times and modes. However switching the plant 
from one mode to another is not free and generates expenditures and, on the other hand, when the 
plant is in a specific mode it provides a profit which depends on that mode. 
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It is well-known that the optimal switching problem is related to systems of reflected backward 
stochastic differential equations (BSDEs for short) with inter-connected obstacles or oblique reflection 
(see e.g. |5J QUI ESI EH E3 EZ]) of the following type: Mi G J := {1, m}, 

' Y s l = + Jf fi(u, r)dr + K i T -K i s - jf Z\dB r , V s < T 

< Y* > max jeJ -i{Y s j - gij (ui, a)}, Vs < T (1.1) 
, Io( Y i ~ max jE j-i{Y s J - gij (u, s))dK l s = 0, 

where J~ % = J — {i}. Actually it is shown in the aforementioned papers that Yq is the optimal profit 
if the plant is in mode i at t = 0, i.e., \/i G J , 

sup E[h 5 (T,u) + f f s (s,u>)ds-J29Z k ,S k+ M + i)l [Tk+1<T ]] (1.2) 



<5=(TVfe,a)fc>oe^ 



fc>0 



where, «4q stands for the set of all admissible strategies starting from mode i at time t = 0, fg the 
instantaneous profit per unit of time when 5 is implemented, <7§ fc ,£ fe+1 is the switching cost from mode 
£fc to mode £k+l (Co = h T o = 0) and finally h s (T) is the terminal profit under 5. 

Additionally, if we are given the solution {Y l )i=±,— ,m then the optimal strategy 5* = (£k> T k) 1S 
uniquely characterized as follows: if we set rj" = 0, and £* = i then for any k > 1 

^1 



r* +1 = inf{t > r*, Y^ = max. £j _ 5 . (Y/ - ^.(u,, t))} A T, 



(1.3) 



G argmax{j, ^ fc = max. 6j . £ . ■■(¥? - ge*j(w,t))}. 
In the case of Markovian setting and when randomness stems from an exogenous Ito diffusion process 



(Xl' x ) s < T , (t,x) G [0,T] x M fe , of the following type: 

dXl' x = b(s, X l s ' x )ds + ct(s, X^dBs, s G [t, T] ; X*'* = x for s < t, 

that is to say, the processes fi(s,w), hi(s,u) and gij(s,co) are deterministic functions of (s , Xl' x (u)) , 
the optimal switching problem is also related to the following system of variational inequalities with 
inter-connected obstacles: V i G J 



min I Vi(t,x) - max (-gij(t,x) + Vj(t,x)), -d t Vi(t,x) - Cvi(t,x) - fi(t,x) 



0; 



(1.4) 



Vi(T,x) = hi(x) 

where C is the infinitesimal generator associated with X l ' x . The process X t,x can be the electricity price 
in the market or the dynamics of factors which determine that price. System fll -4f) is the deterministic 
version of the verification theorem of the optimal switching problem. 

Actually in [15] . the authors have proved that if, mainly, the switching costs satisfy gij(t, x) > 70 > 
then the system (II. 4p has a unique continuous solution (v , . . . , v m ) in viscosity sense and the following 
relationship holds true: 

Vi G J,(t,x) G [0,T] x R k ,Vs G [t,T], = v i (s,X t s ' x ). 
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where {Y l )^j are the processes solution of (jl.ip with /j, and /ij which are deterministic functions 
of (s,Xt' x ). 

So the main objective of this paper is to deal with system (|1,4|) in its general setting, i.e., to study 
the existence and uniqueness of a solution in viscosity sense for the following: V i G J 

min < Vi(t, x) — max (—gij(t, x) + Vj(t, x)), 

-d t Vi(t,x) -Cvi(t,x) -f i (t,x,v 1 (t,x),...,v m (t,x),o T (t,x)D x v%x))} = ; ( L5 ) 
Vi(T,x) = hi(x). 

v 

As previously pointed out, system (jl.5l) is connected with the following system of reflected BSDEs with 
oblique reflection : 

' Y^' x = W{X% X ) + Jf fi(r, Xp x , Y r 1]t ' x , . . . , Y™'^, Z^ x )dr + K^ x - K^ x - jf Z^ x dB r , s < T; 
< Y's ,t,x > max jeJ -i{Y s ht,x - g ij {s,X t s ' x )}, s <T; 

j^{Ys' t,x - max je j-i{Yf' t,3: - g ij {s,X t s ,x ))dKl' t,x = 0. 

(1.6) 

This system has already been studied in [19] in this general form. 

There are at least two motivations for considering this problem. The first one is to extend, as much 
as possible, the Feynman-Kac's representation of the solution of (11 .6f) via the solution of (11. 5p and the 
process X t,x . This issue is very important if we are willing to consider the numerical study of (|1.6I) 
especially in connection with the pricing of gas options in the energy market (one can see e.g. [2] for 
more details). The second one is that solving this problem is a step toward the study of zero-sum 
switching games which are encountered especially in the carbon market by energy firms (see e.g. [2] for 
more details). 

The novelty of this paper lies in the fact that we investigate both existence and uniqueness of a 
continuous viscosity solution of (jl.5p under the following relaxed hypotheses: 

(i) the switching costs are non- negative and satisfy the no free loop condition (see [H3]-(ii) below) 

(ii) for any i £ J , either : 

(a) V k G J~ l , the mappings — > fi(t, x, y 1 , . . . , y k ~ 1 ,y k , y k+1 , . . . , y m , z) are non-decreasing ; 

or 

(b) V k G i7 _i , the mappings y^ — > fi(t, x, y 1 , . . . , y k ~ 1 , y k , y k+1 , . . . , y m , z) are non-increasing. 

In both cases we show existence of a solution for (|l,5p while we have been able to show uniqueness 
only in the case (ii)-(a). 

The closest paper to ours is the one by Elie-Kharroubi [11] where the authors deal also with the 
representation of solution of (II. 6h by the viscosity solution of (jl.5p . However their approach, based 
on the minimal solutions of constrained BSDEs with jumps introduced in [23], is not very satisfactory 
since it induces assumptions on the data of the problem which are either not natural or difficult to 
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verify in practice. Finally note that there are also works related to viscosity solutions of the switching 
problem but their settings and/or approaches are not the same as ours [H El [20j [28] , etc. 

This paper is organized as follows: 

In Section 2 we collect the main assumptions on the data of the problem and we define the notion 
of a viscosity solution for the system (|1.5p . In Section 3, we once more introduce the switching problem 
and provide some results related to solutions of systems of reflected BSDEs with oblique reflection 
which are rather new since they are obtained under weaker conditions than the ones of the literature 
on the subject (see e.g. [19], etc.). In particular, we mainly deal with the non free loop property 
[H3]-(ii) on the switching costs. Those results are basic to deal with the main purpose of this work. In 
Section 4, we provide a comparison result between sub-solutions and super-solutions of the system (11.51) 
in the case when for any i G J the function /j depends on {y % )iej only through y l . We then show, in 
this specific framework of functions /j, that system (II. 5p has a unique continuous solution (v r (t,x))i e j 
which is moreover of polynomial growth. As a by-product, we provide a probabilistic representation 
for the solution of (11.6P via the deterministic continuous functions (v l (t,x))i^j and the process X t,x . 
In Section 5, we deal with the general framework under mainly conditions (i)-(ii) above. Using the 
results of Sections 3 and 4, we construct in each case an approximating scheme which is convergent 
and whose limit is a solution in viscosity sense for system (ll.5p . Finally under condition (ii)-(a), we 
show that system (II. 5p satisfies the comparison property between sub- and super-solutions. Thus under 
conditions (ii)-(a) the solution of system (II. 5|) is unique. ■ 

2 Assumptions and problem formulation 

Let T (resp. k) be a fixed real (resp. integer) positive constant, let J := {1, . . . , m} and let us consider 
the following functions: for i,j G J, 

b : (t,x) G [0,T] x R k ^ b(t,x) G R k ; 
a : (t, x) G [0, T]xl fc 4 a(t, x) G R kxd ; 

fi : (t,x,y\...,y m ,z) G [0, T] x R k+m+d ^ f t (t,x,y l , ...,y m , z) G E ; 
9ij : (t, x) G [0, T]xl fc 4 9ij (t, x)£R(i/ j); 
hi : x G R k ^ hi(x) G R. 

Next let 4> : (t, x) G [0, T] x R k t— > (p(t, x) Elbea function. It is called of polynomial growth if there 
exist two non negative real constants C and 7 such that: 

\<f>(t,x)\ < C(l + |x| 7 ), V(t,x) G [0,T] x R k . 

Throughout this paper, we denote by II 9 the class of functions with polynomial growth and by 
C x ' 2 ([0,T] x R k ) (or simply C 1,2 ) the set of functions defined on [0, T] x R k with values in R which 
are C 1 in t and C 2 in x. □ 
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We now consider the following assumptions: 

[HI]: The functions b and a are jointly continuous and of linear growth in (t, x), and Lipschitz contin- 
uous w.r.t. x, i.e., there exists a constant C > such that for any t G [0, T] and x,x' G M. k 

\b(t,x)\ + \a(t,x)\ <C(l + \x\) and \a(t, x) - a(t, x')\ + \b(t,x) - b(t, x')\ < C\x - x'\. (2.1) 

Throughout this paper we assume that assumption [HI] holds. 

[H2]: for i € J, fi satisfies: 

(i) (t, x) i— > fi(t, x, y 1 , y m , z) is continuous uniformly w.r.t. (j?,z) := (y 1 , y m , z); 

(ii) fi is uniformly Lipschitz continuous with respect to (l/, z) := (y 1 , ...,y m , z), i.e., for some C > 0, 

\fi(t,x,y\...,y m ,z)-f i (t,x,y\...,y m ,z)\ < C^y 1 - ^\ + ■ ■ ■ + \y m - y m \ + \z - z\) ; 

(iii) the mapping (t, x) >->■ x, 0, ... ,0) is fi([0,T] x IR fe )-measurable and of polynomial growth 
i.e. it belongs to U 9 ; 

(iv) Monotonicity : Mi £ J, for any k G J"*, the mapping y fc GR4 x, y 1 , 2/fe+i, 2/m) 
is non-decreasing whenever the other components (t, x, yi, y^-i, 2/fe+i, y m ) are fixed; 

[H3]: (i) gij is jointly continuous in (t,x), non-negative, i.e., gij(t,x) > 0, V(i, x) € [0, T] x R fc and 
belongs to IT 5 ; 

(ii) T/ie non-free loop property: for any (t,x) G [0, T] x IR fc and for any sequence of indices i\,...,ik 
such that i\ = ik and card{i\ : .-.,ik} = k — 1 we have: 

g hi2 (t,x) +g i2i3 (t,x) + ---+g ik _ lik (t,x) + g ikh (t,x) > 0, V(i,x) G [0,T] x M fc . 

As a convention we assume hereafter that gu(t,x) = for any (t,x) G [0, T] x R fc and i € J. 
[H4]: /ij is continuous, belongs to II 9 and satisfies: 

Vx G M, /ij(x) > max (hj(x) — gij((T,x)). ■ 

Next let us introduce the following infinitesimal generator 

Op(t,x) = ^Tr[{a.a T ){t,x)D 2 xx ^{t,x)}+b{t,x) T .D x ^{t,x) (2.2) 

for a function ip which belongs to C 1,2 ([0,T] x R fc ;IR) (Tr is the trace of a symmetric matrix and (.) T 
stands for the transpose). It is associated with a stochastic process which we will describe precisely 
below. ■ 
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In this paper we are concerned with the existence and uniqueness in viscosity sense of the solution 
(v 1 , . . . ,v m ) : (t,x) G [0,T] x R fc \-t (v l (t,x), . . . ,v m (t,x)) G R m of the following system of m partial 
differential equations with inter-connected obstacles: V i G J 

min < Vi(t, x) — max (—gij(t, x) + Vj(t, x)) ; 

-d t Vi(t,x) -Cvi(t,x) - f i (t,x,v 1 (t,x),...,v m (t,x),a T (t,x).D x v i (t,x))} =0; ( 2 - 3 ) 
Vi(T,x) = hi(x). ■ 

To proceed we will precise the notion of a viscosity solution of the system (|2.3[) . It will be done in 
terms of sub- and super-jets. So for any locally bounded function u : (t, x) G [0, T] x R k i— > u(t, x) G R, 
we define its lower semicontinuous (Isc for short) envelope u*, and upper semicontinuous (use for short) 
envelope u* in the following way: 

u*(t,x) = lim u(t , x ) and u*(t, x) = lim u(t,x). 

(t',x')-+(t,x),t'<T (t',x')->(t,x),t'<T 

Definition 1 ; Subjects and superjets 

(i) For a function u : [0, T] x R fc — > R, Isc (resp. use), we denote by J~u(t, x) the parabolic subjet (resp. 

J + u(t, x) the parabolic superjet) of u at (t, x) G [0, T] x R fc , as the set of triples (p, q, M) G R x R k x S k 
satisfying 

u(t', x) > ( resp. <) u(t, x) + p(t' — t)+ < q, x' — x > + 
l<x'-x, M(x' -x)> +o(\t' -t\ + \x' - x\ 2 ) 

where S k is the set of symmetric real matrices of dimension k. 

(ii) For a function u : [0, T] x R fc — > R, Isc (resp. use), we denote by J~u(t,x) the parabolic limiting 
subjet (resp. J + u(t,x) the parabolic limiting superjet) of u at (t,x) G [0,T] x R k , as the set of triples 
(p, q, M) Glxl'x § k such that: 

(p, Q, M ) = lim(f>n, Qn, M n ), (t, x) = lim n (i n , x n ) 

n 

with (p n ,q n ,M n ) G J~u(t n ,x n ) (resp. J + u(t n ,x n )) andu(t,x) =1im n u(t n ,x n ). 

We now give the definition of a viscosity solution for the system of PDE equations with oblique 
reflection fl%3J) . 

Definition 2 : Viscosity solution to Ii2.cl\) 

(i) A function (v\, ...,v m ) : [0, T] x R fe — > R m such that for any i G J , V{ is Isc (resp. use), is called 
a viscosity supersolution (resp. subsolution) to \2. 3\) if for any i G J , for any (t, x) G [0, T) x R k and 
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any (p,q,M) £ J Vi(t,x) (resp. J + v l (t,x)) we have: 

min < V{(t, x) — max (—g%j(t, x) + Vj(t, x)); 

-p - b(t,x) T .q - ±Tr[(aa T )(t,x)M] - fi(t,x, Vl , . . . ,v m ,a(t,x) T .M)} > (resp. < 0); 
Vi(T,x) > (resp. <)hi(x). 

v 

(ii) A locally bounded function (vi, v m ) : [0, T] xI l -> W 11 is called a viscosity solution to 112. 3\) if 
(i>l*> v m *) (resp. (v\ , ...jV^)) is a viscosity supersolution (resp. subsolution) of h2.3\) . ■ 

As pointed out previously we will show that system ()2.3j) has a unique solution in viscosity sense. 
A particular case of this system is the deterministic version of the optimal m-states switching problem 
which is well documented e.g. in [TOj [20] and which we will describe in the next section. 

3 The optimal m-states switching problem 
3.1 Setting of the problem 

Let (£l,J-, P) be a fixed probability space on which is defined a standard d-dimensional Brownian 
motion B = (B t )o<t<T whose natural filtration is (J 7 ® := a{B s ,s < i})o<t<T- Let F = (J : t)o<t<T be 
the completed filtration of (-7"? )o<t<T with the P-null sets of F, hence (Tt)o<t<T satisfies the usual 
conditions, i.e., it is right continuous and complete. Furthermore, let: 

- "P be the <7-algebra on [0, T] x of F-progressively measurable sets ; 

- % 2 ' k be the set of P-measurable, M fc -valued processes w = (wt)t<T such that IE[/ T |w s | 2 cis] < oo ; 

- S 2 be the set of "P-measurable, continuous, M-valued processes w = (wt)t<T such that 
E[sup t<T |u>t| 2 ] < oo. 

The problem of multiple switching can be described through an example as follows. Assume we 
have a plant which produces a commodity, e.g. a power station which produces electricity. Let J be 
the set of all possible activity modes of the production of the commodity. A management strategy of 
the plant consists, on the one hand, of the choice of a sequence of nondecreasing stopping times (r n ) n >i 
(i.e. T n < r n +i and tq = 0) where the manager decides to switch the activity from its current mode to 
another one. On the other hand, it consists of the choice of the mode £ ra , a r.v. T Tn -measurable with 
values in J , to which the production is switched at r n from its current mode. Therefore the admissible 
management strategies of the plant are the pairs (5, £) := ((r n ) n >i, (£n)n>i) for which we require also 
that P[r n < T, Vn > 0] = 0. This set is called of admissible strategies and denoted by V. 

Next, assuming that the production activity is in mode 1 at the initial time t = 0, let (ctt)t<T denote 
the indicator of the production activity's mode at time t £ [0, T], i.e., 

(t) + J2^(r n ,r n+1 ](t). (3.1) 
n>l 
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Finally, for i G J ^ let (ifti(t,uj))t<T be a process of Ti 2 ' 1 which stands for the instantaneous profit 
when the system is in state i and for i, j G J ', i ^ j, let (gij(t,u))t<T be a process of S 2 which denotes 
the switching cost of the production at time t from current mode i to another one j. If the plant is run 
under the admissible strategy (<5, £) = ((r n ) n >i, (£ n )n>l) the expected total profit is given by: 

i-T 



J(6,t) := E[ / i> as (s)ds - Y J 9^-uU T n)^[ Tn< T]}- 
J ° n>l 



Therefore in several works authors are usually interested in either finding an optimal strategy, i.e, a 
strategy (£*,£*) such that J(5*,£*) > J(<5,£) for any (<5,£) G T> ([lOl [20]) or at least in characterizing 
the quantity sup^^g-p J(5,£) in some specific cases (OUS]). This latter quantity is in a way the price 
of the power plant in the energy market. ■ 

3.2 Connection with systems of reflected BSDEs with oblique reflection 

In order to tackle the switching problem described above, we usually relate it to systems of reflected 
BSDEs with oblique reflection which we introduce below in the case we need in order to deal also with 
the system of PDEs (fOD . 

Let (t, x) G [0, T] x W k and (Xt' x ) s <T be the solution of the following stochastic differential equation: 

dX l f = b(s,Xl' x )ds + a{s,Xl' x )dB s , s G [t,T] and X l f = x for s G [0,t]. 

The solution of this equation exists, is unique, since b and a verify [HI], and satisfies: 

Vp > 1, E[sup \XY\ V \ < C(l + |x| p ). (3.2) 

s<T 

Next let us introduce the solution of the system of reflected BSDEs with oblique reflection associated 
with the deterministic functions ((fijiej, {9ij)i,j£j, (hi)iej) introduced in Section 2. The solution 
consists of m triplets of processes ((Y r ' t,x , Z r ' t,x , K r ' t,x ))i(zj, which is denoted, for convenience, by 
((Y\ Z 1 , K l ))i<zj and satisfies: for any i G J, 

Y\ K l G S 2 , Z % G H 2 ' d and K % non-decreasing and Kq = 0; 

Y* = hi(X^) + for, X$\Y r \ . . . , Y™, Zl)dr + K l T - K\ - jj ZidB T , V s < T 
Y* > max jeJ -i{Yi - gyfaX**)}, Vs < T 
{ J T (F; - max jeJ -i{Yi - 9ij (s, xt' x )})dKl = 0. 

We first provide an existence and uniqueness results of the solution of (j3.3|) and some of their 
properties as well. 

Theorem 1 Assume that: 

1) the functions (fi)i£j satisfy (H2)-(ii), (Hi) and (iv) ; 

2) For any i,j G J , the functions (resp. hi) verify (H3) (resp. (H4))- 

Then the system (3. 3]) has a solution ((Y l , Z l , iT l ))i=im- 
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Proof : Since the above assumptions are not exactly the same as the ones of Theorem 3.2 in [19] we give 
its main steps for sake of completeness. 



Step 1: Let us consider the following BSDEs: 

Y eS 2 ,Z e H 2 ' d 

Y s = max i= i im hi(x!f) + /J[max i= i ]m /j](r, Xp x ,Y r , . . . ,Y r , Z r )dr - jj Z r dB r , s <T, 
and 

' Ye5 2 ,Z en 2 ' d 

Y s = minj=i )m hi{X^ x ) + / s T [mini = i im fi](r, Xr' x ,Y r , . . . ,Y r ,Z r )dr - fj Z r dB r , s <T. 



(3.4) 



(3.5) 



Thanks to the result by Pardoux-Peng [23], the solutions of both (|3.4p and ()3.5[) exist. We next introduce 
the following sequences of BSDEs defined recursively by: for any i £ J, Y l '° = Y and for n > 1 and 

s <T, 

' y i,n R i,n e 5 2 ) Z i,n £ y2,d and R i,n non _decreasing; 
Yp n = MX**) + Jf /,(r, Xp*,Y*> n -\ Yr hn '\Yr, Yr +1 ' n ~ 1 , Y^ n -\zT)dr 

+K!f-KT-^Zp n dB r - (3.6) 
Ys ,n > max j€ j- t {Yi' n ~ - g^s, Xi' x )}; 
{ - m^j-^Yr- 1 - gij (s, xt> x )})dKT = 0. 

By an induction argument and the result by El-Karoui et al. (|12j. Theorem 5.2), we claim that the 
processes (Y l ' n , Z l>n , K l,n ) exist for any n > 1. Next using the comparison theorem of solutions of 
BSDEs (see e.g. Theorem 2.2 in [13J ) we deduce that for any i E J, Y 4, ° < Y 1 ' . Now since /j satisfies 
the monotonicity property (H2)-(iv) and using once more the comparison of solutions of reflected BSDEs 
(see e.g. Theorem 4.1 in [12]) we obtain by induction that: 

Vn > and i £ J, Y*' n < Y*' n+1 . 

Since the processes ((Y, Z, 0))j g j- is a solution for the system of obliquely reflected BSDEs associated 
with (([maxpi^/^s,!^,!/ 1 , . . . , y m , z)) ieJ , (max,; =1]m hi(X^, x )) ie j, (gij(s, Xt'*))^^), an induction 
procedure and the repeated use of comparison theorem, which is justified taking into account that fi 
satisfies the monotonicity property (H2)-(iv)), leads to 

Vn > 0, Vie J, Y l ' n < Y. 



Step 2: Using Peng's monotonic limit theorem (see Theorem 2.1 in [25]), we deduce that for any i 6 J , 
there exist: 

(i) a cadlag (for right continuous with left limits) process Y % such that Y t,n /* Y l pointwisely ; 



(ii) a process Z % of T-L 2 ' d such that, at least for a subsequence, (Z* ,n ) n >o converges weakly to Z l in 
/ H 2 ' d and strongly in L p (dt ® cLP) for any p£ [1, 2[ ; 

(iii) a cadlag non decreasing process K l such that for any stopping time r, (Kp n ) n >o converges to 
K* in L p (dP) for any p G [1, 2[. 

Additionally in taking the limit in (13. 6p . the triple of processes (Y 4 , Z*, i<f J ) satisfies: 

Y? = ^(X^) + /i(r, X*' 33 , Y; 1 , . . . , Y™, Z;)dr + itj, - JT* - jf Z\dB r and 

Next let us consider the following m independent reflected BSDEs with a cadlag barrier which, in 
addition, are independent of each other (for existence and uniqueness results for such BSDEs, we refer 
to pi]): 

' Y* = hiixtf) + jf fi(r, X^X, Y;-\ Y r \ Y;+\ . . . , Y r m , Zi)dr + K i T -K i s - jf Z' r dB r , s < T 

Y* > max jeJ -i{Yj - gij (s,xt' x )}, Vs < T 

[ (Y!_ - max i67 _i{y/_ - giJ (s,X^)})dKi = 0, s < T 

(3.8) 

where Y/_ (resp. = lim s/ *4 Y? (resp. ifg), i > 0. Henceforth in using once more comparison 

theorem (e.g. |16j . Theorem 1.5) we have for any i G J , Y l ' n < Y % and then Y % < Y l . On the other 
hand using Ito's formula with ((Y/ - Y/) + ) 2 , t < T, we obtain that Y t { < Y t \ for any t < T. Actually 
this is true because Y l is the minimal solution which satisfies system (13. 7j) . It implies that Y l =Y l , 
and then = Z*, K l = K\ for any i G J". Therefore (Y*, Z\K l ) satisfies not only the system (pTfl) 
but also: 

/ (Y;_ - max {Y/_ - (a, = 0, Vi G J. 

It remains to show that, for each i G J ,Y l is continuous. For this, let us assume that for some i and to, 
AY/ o := (Y/ Q — Yf _) < 0, therefore AiQ > and then from (|3.8[) we deduce the existence of i\ G J~ l 
such that Y£_ = Yj_ - 5iil (t , X*f) and as Y£ > Y^ 1 - 9iil (t , X 1 ^) then AY^ 1 < 0. Repeating this 
reasoning as many times as necessary we deduce the existence of a sequence of indices (i^) such that 
ik 6 J~ ih ~ x and 

Y/^J 1 = Y^_ - g ik _ lik (t ,Xl' o x ). 
As ^7 is finite then there exists a sequence j\, .., j p —i which are all different such that 

g jlj2 (t , X l t ' x ) + ■■■ + g jp _ dl (t , Xl' x ) = 0. 

But this is contradictory with the assumption [H3]-(ii) and then such a to does not exist and Y l is 
continuous for any i G J . Therefore ((Y l , Z % , K l ))i £ j, is a solution for the system of reflected BSDEs 

Q. ■ 
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We now give a remark related to comparison of the solutions of system (|3,3p constructed in Theorem 
[U Its proof is rather easy since an induction argument allows to compare the solutions of the convergence 
schemes. Actually we have: 

Remark 1 : Let (fl)iej (resp^g'^ij^j , resp. (h'^i^j) be functions that satisfy (H2)-(ii), (Hi), (iv) 
(resp. (H3), resp. (H4)) and let {{Y n , Z H , K h ))i £ j be the solution of the system of reflected BSDEs 
associated with {(fl)iej, {g'ij)i,j£j> constructed as in Theorem^ If for any i,j G J we have: 

fi< fl, h < h'i and g i:j > g {j 

then for any % G J , 

In case of uniqueness of the solutions of those systems, this result reduces to the comparison of the 
solutions. ■ 

We now focus on the regularity properties of the solution of system (13, 3j) constructed in Theorem 

m 

Proposition 3.1 Assume the assumptions of Theorem^ are fulfilled. Then there exist Isc deterministic 
functions {v l )i^j, defined on [0, T] x M.-valued and which belong to H 9 , such that: 

Viej,Vse[t,T], y: = v l (s,x t f), 

where {{Y l , Z l , K l ))i^j is the solution of \3. 3\) constructed in Theorem^ 

Proof: Actually under the hypotheses of Theorem [H there exist deterministic continuous with polyno- 
mial growth functions v(t,x) and v(t,x) with values in K such that for any s S [t, T],Y S = v(s,xl ,x ) 
andy s =v{s,X l f) ([13J, Th.4.1). 

Next by induction and thanks to the result by El-Karoui et al. ([E], pp.729), there exist determin- 
istic continuous functions v l ' n (t,x) in the class II s such that for any i G J and n > 0, 

Yi> n = v i > n (s,X t s > x ),Vse[t,T], (3.9) 

the process Y l ' n being defined as the unique solution of (|3SJ) (see step 1, Theoremd]) As Y^ n < Y^ n+1 < 
Y then for fixed i, the sequence (v l ' n ) n >o is non-decreasing and such that v % ' n < v, then it converges 
pointwisely to v l . This latter function is therefore lower semi-continuous on [0,T] x R k , of polynomial 
growth since v < v l < v and finally for any s € [t, T], Y a * = v l (s, X s ' x ). □ 

Remark 2 ; Since v l , i G J , belongs to IT 9 , then classically (see e.g. UEj ) one can show that for any 
i G J , \\Z l \\^2, d (t, x) is also of polynomial growth. ■ 
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We next give a representation result for the solutions of system (|3.3j) and, as a by product, we 
obtain a uniqueness result in some specific cases. Actually let fix it := (ii J ) i=l m in % 2,m and let us 
consider the following system of reflected BSDEs with oblique reflection: 

Vi e J, Vs < T, 



(3.10) 



Y u '\ K u ' i G S 2 , Z u * G V 24 and K"* non-decreasing (K^ = 0); 

= hi(X$*) + jf fi(r, X l r X M, Z^)dr + - KV - jf Z^dB r - 
Ys' 1 > maXj eJ -i{Y?' 3 - g ij {s,X t s ' x )}; 
, Io( YsU ' 1 ~ max igt7 - I {y" J - gij(s,Xt' x )})ds = 0. 

Let s < T be fixed and i G J and let V l s be the following set of admissible strategies : 

V\ := {a = ((r n )„> , (£„)„>o) G 2?,£o = *,t = 0, n > s and £[(A£) 2 ] < oo} 
where for any r G [0,T], A" is the cumulative switching costs up to time r, i.e., 

K ■■=J2 9 ^-i,^{ r n,X t T ^)l [Tn < r] for r < T and A T = KvaA r , Pa.s. 

n>l V 

Therefore and for any admissible strategy a we have: 

n>l 

Note that, by definition of the set D, (r n (uj)) is a stationnary sequence (for almost all oo) and therefore 
the previous sum is finite, P-almost surely. 

Let us now consider a strategy a = ((r n ) n >o, (£n)n>o) G V l s and let (P a ,Q a ) := (P",Qs) s <T be the 
solution of the following BSDE (which is not of standard type) : 

P a is RCLL and £[sup s < T |Pf| 2 ] < oo, Q a G % 2 ' d ; 
^ Pj* = ^ a (X^) + jf / a (r, Q?)dr - jf Q«dB r - {A% - A% V s < T, 

with 

K{x) = ^2h (n (x)l [Tn<T < Tn+l[ and f a (t,x,(C) ieJ ,q) : = ^ / ?n (r, x, (CW, 9)l[r„<r<r„+i[- 

n>0 n>0 

In setting up P a := P a — A a , we easily deduce the existence and uniqueness of the process (P a , Q a ), 
since A a is adapted and 22[(.A!p) 2 ] < oo, and the generator as well as the terminal value of the trans- 
formed BSDE are standard. 

We then have the following representation for the solution of (13.101) which is the main relationship 
between the value function of the optimal switching problem and solutions of systems of reflected 
BSDEs with oblique reflection. This result usually referred as the verification result is not new and has 
been already shown in several contexts and under various assumptions. 
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Theorem 2 Assume that for any i,j G J : 

(i) fi satisfies (H2)-(ii), (Hi) and for any fixed (y, z) the mapping (t, x) h-> fiit, x, y , z) is S([0, T] xR fc )- 
measurable ; 

(H) gij (resp . hi) satisfies (H3) (resp. (H4)). 

Then the solution of system of BSDEs \3.1Q) exists and satisfies: 

Vs < T, V* G J, = ess sup aeVl P«. (3.12) 

77ms i/ie solution of 13.10\) is unique. 

Proof : Thanks to Theorem Q] and considering once more the same assumptions on the functions 
(fi>9ij>hi)i,jeJ the- solution (Y u ' z , Z u '\ K u > l ) ie j of system (j3.10j) exists. Next using that Y u ' 1 is so- 
lution to system (13,101) and following the strategy a G T>\ in (|3.3j) . we obtain: 

Y s u > 1 > h a (X^ x ) + / f a (r, X^, ut, Z?)dr - / Z*dB r - {A% - A") + K$ (3.13) 

</ S J S 

where 

= (7^-7^) if t < r < n, K? = K? n +(K$»-K%) if r n < r < r n+l ; Vr < T, Z r Q = £ Z% n < r<Tn+l[ . 

n>0 

As TiTy > then we have: 

Y?* > If, V a G Pj. 

Note that the right-hand side in f)3. 13j) is not a BSDE, therefore we shall rather consider the equation 
satisfied by Y u ' 1 — P a where the pair (P a ,Q a ) satisfies (|3.1ip . Then using an equivalent change of 
probability we deduce the previous inequality. 

Next let a* = (T*,£*) n >o be the strategy defined recursively as follows: Tq = 0, £g = i and for 
n > 0, 

r* +1 = inf{s > r*, Y^ = max {Y y - g^(s, X$*))} A T 

and 

Let us show that a* G T>\ and, for this, let us first prove that P[t* < T, Vn > 0] = 0. Actually assume 
the contrary i.e. P[t* < T, Vn > 0] > 0. Therefore thanks to definition of r*, we have: 

p Kf, = Y r£ +1 -^ + M + i,Ki^ e; e j-^, vn > i] > o. 

As J is finite then there is a state io G J and a loop io>*i> ■•• 3 *A:,*o of elements of ^7 such that 
cardjio , i\ , ik } = k + 1 and 

p[F «A = y ^ + i _ (r* +1 ,X*f ), Z = 0, . . . , k, (i k+1 = i ), Vn > 1] > 0. 

71+1 71+1 ' 71 + 1 
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Therefore taking the limit w.r.t. n to obtain: 

ppru* = Y ^h +1 _ fl . . +i ( T) x?), l = 0,...,k, (i k+1 = i )\ > 

where r := lim n _ >00 r*. But this implies that 

P[g i0til (r, X?) + ■■■+ g ik , l0 (r, X?) = 0] > 

which contradicts assumption (#3) — (i). Thus we have P[t* < T, Vn > 0] = 0. 

Next it only remains to prove that E[(Aj, ) 2 ] < oo and a* is optimal in T>\ for the switching problem 
(|3.12p . Actually following the strategy a* and since Y u ' 1 solves the reflected BSDE (|3.10p . it yields: for 
any n > 1, 

Y u,i = Y u,C n + r» f^x^^Z? )dr - r ZfdB r - Af* (3.14) 

Js Js 

noting that Kr n — Kr" = holds for any r, r* < r < r* +1 . Taking now the limit w.r.t. n in (|3. 14[) to 
obtain: 

Y s u > % = h a *(X t /)+ f a ,(r,X t r ' x ,^.,Zf)dr- ZfdB r -A$. (3.15) 

J s J s 

But taking into account the assumptions (H4) and (H2)-(ii),(iii) satisfied by hi and fi respectively and 
since it G H 2 ' m and Z a * G V 24 and (Y% G (S 2 ) m , we deduce from (I3T51 that £[(Af ) 2 ] < oo. It 
follows that a* G and Y s u ' 1 = P s a , thus (I3.12|) holds and the solution of (I3.10P is unique. ■ 

Next for it := (it i ) i=1)m G % 2 ' m let us define by 

Hit) := r^!,™ (3.16) 

where if M '*)) i=lm is the solution of system (|3.10j) which exists and is unique if the assump- 

tions of Theorem [2] are fulfilled. Note that when the processes (Y u,l )i = i tm exist they belong to (S 2 ) m 
and then $ is a mapping from T-L 2 ' m to H 2,m . 

The following result, established by Chassagneux et al. [7], shows that $ is a contraction in 1~L 2 > m 
when endowed with an appropriate equivalent norm. Therefore the existence and uniqueness of a 
solution for fj3.3j) is deduced for general functions /, since, contrary to Theoremdl they are not supposed 
to satisfy any monotonicity assumption. Actually we have: 

Theorem 3 Assume that for any i,j G J the following assumptions are fullfiled: 

(i) fi verifies (H2)-(ii),(iii) ; 

(ii) gij (resp. hi) verifies (H3) (resp. (H4))- 
Then we have: 

(a) For any it = (u^l.m, ~& = (^)i=l,m G ^ 2,m , 

rT 

Vi G J, Vs < T, E[|y s ^ - Y^\ 2 } < C(\\lt - lt\\ 2 H 2,m := E[ / ||^ - 4\\ 2 dr\) ; (3.17) 

Jo 
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(b) The mapping $ is a contraction when 7i 2 ' m is endowed with the following equivalent norm: 
||(^)i=l,mlk == M / e* S ( \<\ 2 )ds]}K where it = {u% =l , m € U^ m 

^° i=l,m 

for some appropriate /3q E R. 

Proof : We provide it only for the sake of completeness since it has been already given in [7]. For i € S, 
it and it En 2 > m let us set 

<^(r,^>V) = f i (r,X t r r B ,i$.,z)\/f i (r,Xp*,$,z), r < T, 



and let us consider the solution, denoted by (Y l , Z % , K l )i^j , of the system of obliquely reflected BSDEs 
associated with (<fi(r, Xr' x , z))i e j , {hi)i^j and (gij)i,jeJ which exists and is unique by Theorem[H As 
shown in Theorem O the following representation holds true: 

Vs<T, Yi = ess sn VaeVi P s 

where for a € T>\ the pair of processes (P a ,Q a ) verifies: 

P a is RCLL and £[su Pr? < T \P%\ 2 } < oo, Q a G V 2 ^; 

P% = h a (X f /) + ft Va {r,X^\Q a r )dr - Qf.dB r - (A a T - A"), V?? < T. 

Additionally an optimal strategy a exists i.e. Y£ = P®. Note here that the dependence of P s a on i is 
made through the strategy a which belongs to T>\. Now since for any r <T and z £ W 1 , (fi(r, Xr ,x , z) > 
fi(r, Xr ,x ,u r , z) and (pi(r, Xp x , z) > fa(r, Xp x , v r , z) then by comparison and uniqueness (see Remark 
[T]) we have: 

Y u,i < Y l and r- 1 < Y\ (3.18) 

Next for a £ V*, let (P r a , Q^) r <T be the solution of the non-standard BSDE (fBTTTT) and let (Pf, Q' r a ) r <T 
be the solution of the same non-standard BSDE with generator f a (r, Xr ,x , vt,z). Then we have: 

pa < yu,i < yi = pa p,a < yv,i < yi = pa 

which implies, 

I y s U,i - Y s V,i I < I Ps - P's I + I P 's - P's & I • (3-19) 

But for any rj < T we have: 

pa pa, 

tf{fa(r,Xp x rt,Q°:) V Ur,X^ x y r ,Q- r ) - f a (r,X^ x ,vt r ,Q?)}dr - fi{Q* - Q?}dB r , 
and a similar equation is valid for P3 — P~ a . Next using Ito's formula to obtain: 



\P% ~ P%\ 2 + \Qr ~ Q~r\ 2 dr = -2 fT(Pf - P?){Q~- - Q*}dB r 

+2 ft (pa _ i>«){/ a ( r , Xp*,t£, Q?) V f~ a (r, Xp x ,l$, Qf) - f a (r, X^rfr, Q?)}dr. 
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As for any x,y GM we have \x V y — y\ < \x — y\, then 

\P% ~ P%\ 2 + \Qr ~ Q'r\ 2 dr = -2 f^pa _ pa ){ Q~a _ ga}^ 

+2 \P? - P r d \\h(r, X l r ' x y r , Q«) _ /-( r> X^, ^, Q»)|dr. 
Now classically we obtain the existence of a real constant C > such that: 

l-T 



E[sup |A" - P%\ 2 } < CE[ [ - $fdr]. (3.20) 

V<T Jo 



V<T JO 
}ci r>/a 1 2 



In the same way considering |P" — P^ a | we obtain a similar inequality as (|3.20p where P a is replaced by 
pia Yinathy going back to (|3.19p . squarring and taking the expectation, we obtain the first estimate. 

Let us now show that $ is a contraction. Let (3 > and let us make use of Ito's formula to obtain: 

e ^|pa _ pa|2 + JT & p r \Qa _ Q a^ dr = _ 2 jT (pa _ pa)(ga _ Q§ dBr _ £ ^ e ^|pa _ pa^ 

+2 e^(P r s - if){/ Q (r, X? 8 , «», Q«) V / 5 (r, X*' x , u», Q*) - / g (r, X l r ' x , 5», Q?)}dr, r/ < T. 

But 

|/ a (r, Q«) V / a ( r , X^, Q*) _ / a ( r , X^, ^, Q«| < C {\^. - 4\ + |Q« - Q*|) 

where C is the Lipschitz constant of /. Therefore taking expectation in the previous equation and using 
both inequalities 2Cxy < (Cx) 2 + y 2 and 2xy < + \ffiy 2 for any x, y G M, we obtain: V77 < T, 

E[e^|P a - P^| 2 ] < (C 2 + Cy^ - /3)E[ / T e^|^? - P r a | 2 dr] + -^=E[ f 

Jt] VP Jri 

But the same estimate can be obtained for E[e^|P| - P^\ 2 ]. Taking p > C 2 + Cyfp and going back 
to ([3TT31) to obtain: 

2C f' T 

E[e^ s \Y s u ^ - Y^\ 2 } < -^E[ / ||^ - ^ r \ 2 dr}. (3.21) 
VP Jo 

Next summing for i = l,m in fj3.21f) and integrating w.r.t. dt we obtain that: 



t-T 

\vJr — vt\ 2 dr] 



I— > — >i|2. 



|u r — v r II dr]. 



- *C^)II? := E[ / E " Y s' l \ 2 )ds] < -TS-n \ 

J ° i=l,m VP JO 

Choosing now j3 = (3q > 2 max{(2CT) 2 , C 2 + Cyf]3} yields that $ is a contraction in the Banach 
space (% 2,m , ||.||/3 ), therefore it has a fixed point (l rl )i=i,m which can be chosen continuous since 
<3?((Y*)j = i jm ) 6 (5 2 ) m . Thus the system of reflected BSDEs with interconnected obstacles has a unique 
solution. ■ 

Remark 3 Let {Y ifi ) i(L j be fixed processes ofH 2 > m and for n > 1 let us set (Y i ' n ) ieJ = $((Y i ' ft - 1 ) i6<7 ). 
T/ien the sequence ((Y l ' n )i^j) n >Q converges in (J-L 2,m , ||.||) to i/ie unique solution of the system of 
reflected BSDEs associated with {{fi)i^j,(gij)i,jej,{hi)i & j) since $ is a contraction in (H 2,m , ||-||s Q ) 
and the norms ||.||/3 and \\.\\ are equivalent. ■ 



16 



4 Uniqueness of the solution of the system of PDEs 



In this section we deal with the issue of uniqueness of the solution of system (12. 3D and to do so, we first 
establish an auxiliary result which is a classical one in viscosity literature (see e.g. [26], pp. 76). 

Lemma 4.1 Let (vi(t,x))i = \ m be a supersolution of the system \2. 3\) . then for any 7 > there exists 
Ao > which does not depend on 9 such that for any A > Ao and 9 > 0, the m-uplet (vi(t,x) + 
9e~ xt \x\ 2 ' y+2 )i = i jm is a supersolution for h2.S\ ). 

Proof: Without loss of generality we assume that the functions vi,... ,v m are Isc. For sake of conve- 
nience, we do not use the previous definition of a supersolution but an equivalent one (see e.g. [8]). Let 
i £ J be fixed and let tp G C 1,2 be such that the function tp — {v l + 9e~ xt \x\ 2 " t+2 ) has a local maximum 
in (t,x) which is equal to 0. As (i>i)i=i, m is a supersolution for (|2.3p . then we have: Vi G J, 



min < Vi(t,x) — max (— gij(t, x) +Vj(t,x)) ; 

-8 t {ip{t,x) - 9e- xt \x\ 2l+2 ) ~ \Tria.o- 1 ' (t,x)D 2 x (i P ( y t,x) - 9 e - xt \x\ 2l+2 )) 
-b(t,x) T .D x (ip(t,x) - 9e- M \x\ 2 ^+ 2 ) - fi(t,x, {v 1 ,.. . , v m ){t, x), a T (t, x)D x (tp(t, x) - 9e~ xt \x\ 2 ^+ 2 ))} > 



which implies that 



>i(t,x) + 9e~ xt \x\ 2 ^ +2 ) - max (- gij {t,x) + {vAt, x) + 0e _At |x| 2 T +2 )) 



Vi(t,x) — max (—gij(t,x) + vAt,x)) > 0. 



(4.1) 



On the other hand: 



-d t (ip(t,x) - 9e- xt \x\ 2l+2 ) ~ \Tr{o-.a T {t,x)D 2 xx (<p{t,x) - 9 e - xt \x\ 2 ^ +2 )) 

-D x (ip(t,x) -0e- xt \x\ 2 f +2 ).b{t,x) - f l (t,x,(v l ,...,v m ){t,x),o- T (t,x)D x (cp(t,x) - 9e- xt \x\ 2 "i +2 )) > 
and then 

—d t (f{t,x) - \Tr(a.a T (t,x)D xx (p(t,x)) - b(t,x) T .D x tp(t,x) 

-fi(t, x, (v% x) + 9e- xt \x\ 2 ^ +2 ) l=1 , m , a T (t, x)D x <p(t, x)) 
> 9Xe- xt \x\ 2 ^ +2 - \9e- xt Tr{o.o T {t,x)D 2 xx \x\ 2 ~< +2 ) - 9 e - xt D x {\x\ 2 ~> +2 ).b{t,x) (4.2) 
+[fi(t, x,(v\..., v m )(t, x),a T (t, x)D x (<p(t, x) - 9e~ xt \x\ 2l+2 )) 

-fi(t, x,(v\..., v m )(t, x),a T (t, x)D x tp(t, x))}. 

But the last term in the right-hand side of this latter inequality is equal to 

9e- xt Cl xAX .a T (t,x)D x (\x\ 2 ^ 2 ), 

where CI e A is bounded by a constant independent of 9 since the function /j is uniformly Lipschitz 
w.r.t. z. Therefore, taking into account the growth conditions on b and a, there exists a constant 
Ao £ K + which does not depend on 9 such that if A > Ao, the right-hand side of (|4.2j) is non-negative. 
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Henceforth, noting that i is arbitrary in J together with (|4.ip . we obtain that (v % + 9e A '|x| 27+2 )j = i jm 
is a viscosity supersolution for (12. 3|) , ■ 

We now establish the comparison property between supersolutions and subsolutions of (|2.3p in the 
case when fa does not depend on (y 1 , . . . , y* -1 , . . . , y m ) for any i G J . Actually let us introduce 
the following assumption on the functions /j's. 

[H5]: For any i G J, the function fi does not depend on (y l , . . . ,y l ~ 1 ,y %+1 , ■ ■ ■ ,y m ). Note that this 
assumption replaces assumption (H2)(iv): this last one does not make sense any more when fi depends 
only on y l , z l . 

Proposition 4.1 Assume both (H3) and (E~4) and let suppose that the functions fi, i G J , verify (H2)- 
and (H5). Let (u l (t, x))i=i >m (resp. (v l (t, x))i=i tm ) be a subsolution (resp. a supersolution) of 
the system 12. 3\) which belongs to II 9 , then for any i G J , we have: 

V(t,x) e [0,T] x R fc ,^(f,a;) < v%x). 

Proof: First w.l.o.g we assume that Ui (resp. w l ) is use (resp. Isc) for any i € J . Next let 7 > and 
C be such that that for any i £ J we have: 

\ui(t,x)\ + a?) I < C(l + |x| 7 ), V(t,x) G [0,T] x R k . 

For sake of clarity, the proof is divided into two steps. 

Step 1: To begin with we additionally assume that the functions fi, i G satisfy: 

3A < s.t. Vi, x,z, \/ u > v, fi(t, x, u, z) — fi(t, x, v, z) < \(u — v). (4-3) 

According to the previous lemma we know that for any 9 > and A large enough (vi(t, x) + 
6e~ xt \x\ 2 ' yJr2 )i = i jrn is also a supersolution for (12. 3j) . Therefore it is enough to show that for any i G J\ 
we have: 

V(t, x) G [0, T] x R\ u% x) < v l (t, x) + 6e- xt \x\ w , 

since in taking the limit as 6 — > we obtain the desired result. So let us set w t,e,x (t, x) = v l (t,x) + 
#e~ A *|x| 27+2 , (t,x) G [0, T] x M fc and we still denote w l ' e,x by w*. Next assume there exists a point 
(i, x) G [0, T] x M fc such that for i £ J: maxjgj- (u* (t , x) — w l (t, x)) > 0. Next using the growth condition 
there exists R > such that: 

V(t,s) G [0,T] x M fc s.t. |x| > i?, - ^(t,x) < 0. 
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Taking into account the values of the subsolution and the supersolution at T, it implies that 

< max (M;)g[0)T]xRfc maxiej^i^x) - ta'ftx)) = 

max (i,i)6[o,r[xB(o,fi) max i&J (u l (t,x) -w l (t,x)) = max ie:7 (u J (r, x*) -w l (t*,x*)), 

where B(0,R) is the open ball in R k centered in and of radius R and (t*,x*) G [0, T[xB(Q, R). 
Now let us define J as: 

J := {j G J,Uj(t*,x*) -w 3 \t\x*) = max(u k (t*,x*) - w k (t* ,x*))}. (4.5) 

First note that J is not empty. Next for j G J and n > 1, let us define: 

& n {t,x,y) :=u j (t,x)-w^t,y)-<p n (t,x,y), (t,x,y) G [0,T] xl a (4.6) 

where: cp n (t,x,y) := n\x - y\ 2 ~i +2 + \x - x*\ 2 + |t-t*| 2 . Now let (t n ,x n ,y n ) G [0,T] x ^'(O,^) 2 be such 
that 

& n (tn, x n , y n ) = ™ax &> n (t, x, y), 

(t,x,y)£[0,T]x B' (0,R) 2 

which exists since $n is use (B' (0, R) is the closure of B(0, R))- Then we have: 

& n (t*,x*,x*) = Uj(t*,x*)-w j (t*,x*) < Uj(t*,x*)-w j (t*,x*) + ip n (t n ,x n ,y n ) < Uj(t n , x n ) - w j (t n ,y n ). 

(4.7) 

The definition of ip n together with the growth condition of Uj and w 3 implies that (x n — y n ) n >i converges 
to 0. Next for any subsequence ((t ni , x ni ,y ni ))i which converges to (t, x,x) we deduce from (j4.7j) that 

Uj(t* , X*) — W 3 (t* , X* ) < Uj(t, x) — w 3 (t, x), 

since Uj is use and w 3 is Isc. As the maximum of Uj — w 3 on [0, T] x B'(0,R) is reached in (t*,x*) 
then this last inequality is actually an equality. It implies, from the definition of <p n and (|4.7j) . that the 
sequence ((J>m%niyn))n converges to (t*,x*,x *) from which we deduce 

n\x n - ynl 2 ^ 2 ->-n and (uj(t n , x n ),w 3 (t n ,y n )) ->„ (uj(t* ,x*),w 3 (t* ,y*)). 

Actually this latter convergence holds since from (|4.7p we first obtain, 

Uj(t*,x*) - w j (t*,x*) < )imuj(t n ,x n ) - lim w 3 (t n ,y n ), 

n n 

whereas the fact that Uj (resp. Wj) is use (resp. Isc) gives 

limu,(t n ,x n ) - lim-w^tn,^) < Uj(t*,x*) - w 3 (t*,x*). 
n n 

All these inequalities imply that 

lim Uj(t n , x„.) = limuj(t n , x n ) = Uj(t* ,x*) and lim w 3 (t n , x n ) = limw 3 (t n , x n ) = w 3 (t* , x* ) . 

n n n n 
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Next as in [21] . let us show by contradiction that for some k G J we have: 

u k (t*,x*) > max (uj(t*,x*) - g kj (t* , x*)). 
j&J k 

Actually suppose that for any k G J we have: 

u k (t*,x*) < max (uj(t*,x*) - g kj (t* , x*)), 
then there exists j G ,y~ fc such that 

u k {t*,x*) - v,j(t*,x*) < -g kj (t*,x*). 
But w k is a supersolution of (|2.3p . therefore we have 

w k (f,x*) > W J(t*,x*)-g kj (t*,x*) 

and then 

u k (t*,x*) - Uj(f,x*) < -g kj (t*,x*) < w k (t*,x*) - w j (t*,x*). 

Therefore 

u k (t* , x* ) - w k (f , X* ) = Uj (t* ,x*)- W j (t* , X* ) 
which implies that j also belongs to J and 

u k (t*,x*) - Uj(t*,x*) = -g kj (t*,x*). 

Repeating this procedure as many times as necessary and since J is finite we get the existence of a 
loop of indices i±, ...,i p ,i p+ i of J such that i\ = i p+ \ and 

gi^iiit*, x*) H h gi p ,i p+1 (t*,x*) = 0. 

But this contradicts the assumption (H3) on gij, i,j G J, whence the desired result. □ 
To proceed let us consider k G J such that: 

u k (f,x*) > max ( Uj (t*,x*) - g kj (t*,x*)). (4.8) 

As the functions Uj, j G are use and are continuous, then there exists p > such that for 
(t,x) G B((t*,x*),p) we have u k (t,x) > maxj £ j-k(uj(t, x) — g k j(t,x)). Next and by construction it 
holds that (t n , x n , u k (t n , x n )) n — > n (t* , x* , u k (t* , x*)) and once more since Uj is use then for n large 
enough we have: 

u k (t n ,x n ) > max (uj(t 

ni%n) 9k j 

(t n ,x n )). (4.9) 
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Now applying Crandall-Ishii-Lions's Lemma (see e.g. [8] or p3], pp.216) with <& n (note that k £ J 
and ($JB) is satisfied) in (t n ,x n ,y n ), there exist (p™,q™,M™) G J 2 > + u k (t n ,x n ) and (K,,CM™) G 
J 2 ~w k {t n: y n ) such that: 

PS -Pw = d t pn(t n ,x n ,y n ), q% (resp. g") = d x ip n (t n , x n , y n ) (resp. = -d y ip n (t n ,x n ,y n )) and 

M™ \ 1 9 



o -m 



< A n + — Ai (4.10) 



where ^ n = D 2 x ^ip n (t n , x n ,y n ). But 

d t (p n k(t,x,y) = 2(t - t*), d x ip n (t,x,y) = 2(7 + l)n(x - y)\x - y| 27 + 2(x - x*) and 
d y ip n (t,x,y) = 2n(7 + l)(x - y)\x - y\ 21 . 

On the other hand, 

I Dl x tp(t n ,x n ,y n ) Dl y <p(t,x,y) 



with 



DL(f(t,x,y) D 2 (p(t,x,y) 



D 2 xx <p(t, x, y) = 2n( 7 + l)\x - y\ 2 ^ 2 {\x - y\ 2 I k + 2 7 (x - y)(x - y) T } + 2I k , 
D 2 yy v(t, x, y) = 2n( 7 + l)\x - y\ 2 ^ 2 {\x - y\ 2 I k + 2 7 (x - y)(x - y) T } 

and finally 

D 2 xy <p(t, x, y) = -2n( 7 + l)|x - y\ 2 ^ I k - 4n 7 ( 7 + l)(x - y)(x - y) T \x - y\ 2 ^ 2 . 

As {uiji^j (resp. {w l )i & j) is a subsolution (resp. supersolution) of (|3.3p and taking into account (|4.9|) . 
we obtain: 

~Pu - b(t n ,x n ) T .q™ - ^Tr[((Tcr T )(t n ,x n )M™] - fj(t n ,x n ,Uj(t n ,x n ),cr(t n ,x n ) T .q™) < 

and 

-pl-b(tn,y n ) T -qZ ~ ^Tr[{aa T )(t n ,y n )M™} - f j (t n ,y n ,w j (t n ,y n ),a(t n ,y n ) T .qZ) > 0. 
Making the difference between those two inequalities yields: 

-(PS -Pi) ~ (b(t n ,x n ) T .qZ - b(t n ,y n ) T .<£) - ±Tr[{aa T (t n , x„)M« - cr^^, y n )M£}] 
-(fj(t n , X n , Uj(t n , x n ), a J (t n , x n ).q£) ~ fj(t n , Vn, w 3 (tn, Vn), <? T (t n , y^-Qw)) < °- 

But 

lim {\b(t n ,x n ) T .qZ - b(t n ,y n ) T .qZ)\ + \a(t n , x n ) T .q% - a(t n ,y n ) T .<£)|} = 0. (4.12) 

n— >oo 

As usual taking into account (I4.10P we have: 

h^Tr[{aa T (t n ,x n )M: - aa T (t n ,y n )M^}} < 0. (4.13) 
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Finally taking the limit in (|4.11j) and using the assumption (|4.3p to obtain: 

-\( Uj (t*,x*) -w j {t*,x*)) < 0, 
which is contradictory with (|4.4p and then for any j £ J we have Uj < w J . ■ 
Step 2: The general case. 

Once more let (uj)j^j (resp. (w^j^j) be a subsolution (resp. supersolution) of (|2.3j) . For j E J let us 
set Uj(t,x) = e xt uj(t,x) and w j (t,x) = e xt w j (t,x), (t,x) € [0,T] x M fc . Then {uj)jej (resp. (w j )j e j) 
is a subsolution (resp. supersolution) of the following system of variational inequalities with oblique 
reflection: for any i G J, 

min < Vi (t, x) — max (— e xt gij(t, x) + Vj(t, x)) , 
{ jeJ-* 

-d t Vi(t, x) - Cvi(t, x) + \vi(t, x) - e xt fi(t, x, e^v^t, x), e - xt a T (t, x).D x v% x))} = 0; 

Vi(T,x) = e XT hi(x). 

(4.14) 

Actually let i E J and let ip(t, x) be a C 1,2 -function such that ip — Ui has a minimum at (t, x) and 
yj(t, x) = Ui(t,x). Therefore e~ xt Lp — m has a minimum at (t,x) and e _A V(*> x ) = Ui(t,x). As Uj is a 
subsolution then Ui(T,x) < e XT hi(x) and, 

min < Ui(t, x) — max {—Qij{t, x) + Uj(t, x)) , 

-d t (e- xt ip(t,x)) - C(e~ xt ip)(t,x) - fi(t,x,e- xt tp(t,x),a T (t,x).D x e- xt <p(t,x))} < 0. 

(4.15) 

Now if 

Ui(t, x) — max (—gij(t, x) + Uj(t, x)) < 0, 

then we have 

Ui(t,x) — max (— e t gij(t, x) + Uj(t,x)), 

and the viscosity subsolution property is satisfied. If not, i.e., U{(t, x) — max (—gij(t, x) + uj(t, x)) > 
then: 

-d t {e- xt y{t,x)) -£(e- xt <p)(t,x) - fi(t,x,e- xt <p{t,x),a T (t,x).D x e- xt <p(t,x)) < 
which implies 

-d t (p(t,x) + \(p(t,x) - £<p(t,x) - e xt fi(t,x,e~ xt Ui(t,x),e~ xt a T (t,x).D x (p(t,x)) < 0. 

Therefore once more the viscosity subsolution property is satisfied. As i is arbitrary in J , then (ui)i£j 
is a viscosity subsolution for (14.15p . 

In the same way one can show that (w J )j£j is a viscosity supersolution of (I4.15p . whence the claim. 
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Next foiiGj let us set: 

Fi(t, x, u, z) = -Xu + e xt fi(t, x, e~ xt u, e~ xt z). 

Taking now A = (1 + maxj = i m C l ), where C l is the Lipschitz constant of /j w.r.t. to u, to obtain 
that for u>v, Fi(t,x,u, z) — Fi(t,x,v, z) < —(u — v). It means that Fi satisfies the assumption (|4.3p . 
Therefore and according to the result proved in Step 1, for any j G J } we have Uj < w 3 and also 
Uj < w J . The proof of the proposition is now complete. ■ 



Next, thanks to Proposition 14,11 we classically deduce both uniqueness and continuity results of any 
solution of (12. 3ft which belongs to n 9 . Actually if (uj)j 6 j- is a solution then (u*)i^j (resp. (uj*)jej r ) is 
a subsolution (resp. supersolution) for the system (|2.3h in the class II 9 , then we deduce that u* < 
and then u* = = Ui, for any i £ J . Whence the continuity of (ui)i e j. To sum up, we have: 

Theorem 4 Assume that Assumptions (H2)-(i),(ii) and (H5) for fi, i £ J , are fulfilled. Then: 

(i) The system of variational inequalities with inter- connected obstacles \2.3\) has at most one solu- 
tion in the class H 9 ; 

(ii) If the solution in H 9 exists, it is necessarily continuous. ■ 

5 Existence of a solution for the system 

5.1 Case 1: fi depends only on (y\z l ). 

In this specific case we have the following existence result: 

Theorem 5 Under (H2), (H3), (H4) and (H5), the following system of variational inequalities with 
inter- connected obstacles 

min < Vi(t, x) — max (—gij(t, x) + Vj(t, x)) , 

-d t Vi(t,x)-£ Vi (t,x)- fi^x^Viit^.o-^it^.D^it^x))} =0; ( 5>1 ) 
Vi(T,x) = hiix) 

v 

has a unique continuous solution (v l )iej in the class II 9 . 

Proof: First note that the hypothesis (H2)-(iv) does not make any sense in consideration with (H5). 
Now let (v l )i£j be the functions constructed in Prop. 13.11 which are associated with the solution of the 
system of reflected BSDEs with inter-connected obstacles associated with {{fijiej-, (hi)iej> (9ij)i,jej), 
which both exist under (H2)-(H5). The functions v l , i G J , are of polynomial growth, thus locally 
bounded. Next let us show that they are viscosity solutions for the system (|5.ip . 
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For any i G J', v l is Isc, then v % = v\. So let us show that the m-uplet (v l )i G j is a viscosity 
supersolution to (I5.ip . First note that for any i S J, 



u l = lim /• v l,n , 

n— >oo 



(5.2) 



where v l ' n , for n > 1, is defined in (|3.9|) . By El-Karoui et al.'s result ([12], Thm. 8.5), v l,n is a viscosity 
solution of the following variational inequality or PDE with obstacle: 

mm{v l ' n (t, x) — max J - g jr-i[ t, "'' n ~ 1 (^ a; ) — gij(t,x)]; 

-d t v^ n - Cv^ n {t, x) - fi(t, x, v^ n (t, x),a{t, x) T J D x v i ' n (t, x)} = 0; (5.3) 
{ v i < n (T,x) = h i (T,x). 

Now let us fix i £ J, let (t,x) £ [0,T[xR k and (p,q,M) £ J-?/(i,x). By $x2$) and Lemma 6.1 in [8], 
there exist sequences 

rij -> oo , (pj , , Mj ) £ J" ■u*' nj (i,- , Xj ) 

such that: 

{tj,Xj,v hn i(tj,Xj) y pj,qj y Mj) ->• (t,x,v l (t,x),p,q,M). 
Now from the viscosity supersolution property for v l ' n i we have: 

-Pj - b(tj,Xj) T qj - irr(cro- T (tj,Xj)Mj) - fi(tj ,Xj,v l ' n i (t,x),a(tj ,Xj) T qj) > 0, 

and taking then the limit as j — > oo we obtain: 

-p - b(t, x) T q - ^Tr(aa T {t, x)M) - fi(t, x, v% x),a(t, x) T q) > 0. 

As v % (t,x) > maxjfzj-i^ (t, x) — gij(t,x)) and v i (T,x) = hi(x) then v l is a viscosity supersolution for 
the following PDE with obstacle: 

mm{v' l (t, x) — max je j- 1 [v 3 (t, x) — gij(t,x)]; 

-dtv* - Cv l {t, x) - hit, x, v% x),a(t, x^D^it, x)} = 
[ v i (T,x) = h i (x). 

Finally as i is arbitrary in J then the m-uplet (y , . . . ,v m ) is a viscosity supersolution for the system 
of variational inequalities (|5.ip . □ 

Next let us show that (v l *)i^j is a subsolution for (|5.ip . First let us show that for any i G J, 
v % *(T,x) = hi(x). To begin with we are going to show that: 



Actually 



mm{v l *(T,x) - hAxV^iT^x) - max (v j *{T, x) - gij (T, x))} = 0. 



v i *(T,x)= lim vHt'x')> lim v^it' x'), for any n > 

{t',x')-MT,x),t'<T (t',x')-+(T,x),t'<T 
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therefore 

> v ^ n (T,x) = hi(x) (5.4) 
since v %,n is continuous and at t = T it equals to /i,(x). On the other hand for any (t,x) we have: 

?/(i,x) > max (v j (t, x) - gij(t,x))}, 



then 



which with ()5.5[) imply that: 



v*(T,x)> max(vi*{T,x)- gij (T,x)), (5.5) 



mm{v l *{T,x) - hdx); v l *(T, x) - max (v^*(T,x) - gij {T,x))} > 0. (5.6) 



Let us now show that the left-hand side of (|5.6p cannot be positive. We first follow the same idea 
as in [3j. So let us suppose that for some xq, there is e > such that: 

mm{v l *(T,x ) - hi(x );v l *(T,x ) - max (v j *(T,x ) -gij(T,x ))} = 2e, 

and let us construct a contradiction. Let (t k , x k ) k >\ be a sequence in [0, T] x M fc such that: 

(t k ,x k ) -)• (T,x ) and v l (t k ,x k ) -)• t>"(T, x ) as /c ->• oo. 



Since v 4 '* is use and of polynomial growth and taking into account of (|5.2p . we can find a sequence 
(£ n ) n > of functions of C 1,2 ([0, T] x R fc ) such that p n -> v^* and, on some neighbourhood B n of (T, xq) 
we have: 

mm{g n (t,x) - hi(x),Q n (t,x) - max (v j *(t,x) - gij(t,x))} > e, V(x,t) e B n . (5.7) 

After possibly passing to a subsequence of (t k ,x k ) k >i we can then assume that it holds on B^ := 
[£fc,T] x B(x k ,5n) for some <5jJ G (0,1) small enough in such a way that C B. Now since t>*'* is 
locally bounded then there exists £ > such that < £ on B n . We can then assume that g n > — 2£ 
on B„. Next let us define: 



4C|x-x fc |2 

Note that > g n and 



&(t, x) := Q n (t, X) + ^ {6n y + VT~t. 



- x) < -C for (t,x) G [i fc ,T] x SB(x fcj <^). (5.8) 

Next since dt(y/T — t) — > — oo as £ — > T, we can choose t k large enough in front of <5j? and the derivatives 
of g n to ensure that 

-£g(t,i)>0on4 (5.9) 
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Next let us consider the following stopping time 6^ := inf{s > t k , (s, X t s k ' Xk ) G B^ C }AT where B^ c is the 
complement of B%, and # k := inf{s > t k ,v i (s,X t s k > Xk ) = max jeJ - t (v j {s, Xt k ' Xk ) - gij (s, X l s k ' Xk ))} A T. 
Applying now Ito's formula to the process X s )) stopped at time 9^ A d k and taking into account 

(pT7|) . ([55D and ([53D to obtain: 

ejf(t fc ,x fc ) = mM A ^,^|' A 1 fc ) " f: A " k Cft(rX h ' Xk )dr] 

> E[{(^*(^,xJ|^) + 0%<t] + (e + ^(r,x*f 9, *))i [< , J , =21 }i [fl *< tfk] 

> E[{(^(0*,X*|«) + C)1 K <T] + (e + ^(T,!^^))!,^^}!^^^, 

+{ e + max ie ^(^(^,X*f fc ))}lp fc «,*]] since > v* 

>E[v*(e k n A$ k ,XlV Xk k )] + (/\e 

= E[«*(t fc , a; fc ) - jJ A,?fc X**'**, X s tfc '**), + C A e 

since on [ifc,i?fc], dK l,t ' x = 0. Finally since and ||Z 4 ||-^2,d(t, x) belong to li 9 (see Prop J3TT1 and Remark 
[2]) and taking into account (|3.2p and assumption (H2)-(iii), we easily deduce that 

lim E[ / fi(s, X^'^ ^'(s, Xl k ' Xk ), Z l / k ' Xk )ds} =0. 
k ^°° Jt k 

Therefore taking the limit in the previous inequalities yields: 



lim Q%(t k ,x k ) = lim g n (t k ,x k ) + ^T-t k = Q n {T : x ) > lim v\t k , x k ) + C A e = v H (T, x ) + C A e. 

k— >oo k— >oo k — ^oo 

But this is a contradiction since g n — > v l * pointwisely as n — > oo. Thus for any (t,x) G [0, T] x M fc we 
have: 

min{?/*(T,x) - ^(x);^*(r,x) - max (v j *(T,x) - gi AT,x))\ = 0. 

Let us now show that v l *(T,x) = h{{x). So suppose that v l *(T, x) > hi(x), then for for some 
j E J~ % we have 

v i *(T,x)=vi*(T,x)-g ij (T,x). 
But once more we have v J *(T,x) > hj(T,x). Otherwise, i.e. if v J *(T,x) = hj(T,x), we would have: 

hi(x) < hj(x) - gij(T,x), 

which is contradictory with (H4). Therefore there exists I £ J~ 3 such that: 

v j *(T, x) = v*(T, x) - gji(T, x) and then v l *(T, x) = v*(T, x) - gjt(T, x) - gij(T, x). 

Repeating this reasoning as many times as necessary we obtain a sequence of different indices i\, 
such that 

v ll *(T,x) = v ll *(T,x) - (g hi2 (T,x) H h g k _ xil (T, x) + g^T, x)), 
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which is contradictory with (H3)-(ii). Thus for any i G J we have: 

Vx G R k , v^(T,x) = hi{x). □ 

Let us now show that (v l *)i^j is a subsolution to (|5.ip . First note that since v l ' n /* v l and v l ' n is 
continuous then we have (see e.g. [26], pp.91) 

= lim supV' n (t,x) = lim sup v^ n (t',x'). (5.10) 

Next let us fix % G J and let (t,ar) G [0,T[xK fc be such that 

?/*(t,x)- max (v £ *(t,x) - gu(t,x)) >0. (5.11) 
Let (p,q,M) G J + ?/*(i,x). By (|5.10p and Lemma 6.1 in JS], there exist sequences 

rij -)■ oo, (pj, qj, Mj) G J^v 1 '^ (tj,Xj) 

such that: 

lim (i,- , xj , u*' 7 ^ (£,• ,Xj), pj ,qj,Mj) = (t,x, v l * (t,x), p, q, M) . 
Now from the viscosity subsolution property for v l ' nj at (tj, Xj) (see I5.3p . for any j > 0, we have: 
mm{v % > n J (tj,xj) -max feJ - ! (/ ,l r 1 (t J)2 ; j ) ~ 9ie(tj,Xj)); 

(5.12) 

Next the definition of t/* implies that 

v e *(t,x) > lim sup v i,n >(tj,Xj), 

j^oo 

therefore by (|5.1ip . there exists jo > 0, such that if j > jo we have 

v^i (tj ,Xj)> max (t/» n i (tj , xj ) - Qij (tj ,xj)). 

Then ()5.12p implies that, for any j > jo, 

-Pj - b(tj, Xj ) T qj - ^Tr(aa T (tj, X j)Mj) - fi(tj,Xj,v l ' n J (tj,Xj),a(tj,Xj) T qj ) < 0. 

Taking the limit as j — > co we deduce that 

-p-b(t,x) T q-^Tr(aa T (t,x)M) - f i (t,x,v i *(t,x),a(t,x) T q) < 

since /j is uniformly continuous in (t,x) and Lipschitz in (y l ,z). Then 

mm{v r *(t,x) -max ieJ -i(v e *(t,x) - gu(t,x)); 

-p-b(t,x) T q - \Tr(aa T (t,x)M) - f i (t,x,v i *(t,x),a(t,x) T q)} < 

which means that v l * is a viscosity subsolution for (|5.ip . Thus the m-uplet (v l )i^j is a solution for 
(|5.ip and Theorem |4] implies that it is continuous and unique. ■ 
As a by-product we obtain: 
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Corollary 1 Under (H1)-(H5), there exist deterministic continuous functions (v l (t,x))i e j which be- 
long to IP unique solution of 15\ ) and such that the unique solution of the system of reflected BSDEs 
with inter- connected obstacles associated with ((ft)iejj(hi)ieji(9ij)i,jej) has the following representa- 
tion: 

Vi G J,V(t,x) G [0,T] x R fc ,Vs G [t,T], Y l s ' l ' x = ?/(s, A^' x ). ■ (5.13) 
5.2 Case 2: Existence and uniqueness in the general setting 

In this section we deal with the issue of existence and uniqueness of the solution of (|2,3p in its general 
form, i.e., when the functions fi, i G J , depend not only on y l but also on the other components 
(yi, . . . , y m ). More precisely, we successively consider the cases when either fi or — fi satisfies (H2)-(iv) 
for any i G J . Considering the first case, we have: 

Theorem 6 Under (H2), (H3), (H4), the system of variational inequalities with inter- connected ob- 
stacles \2. 3\) has a continuous solution (v , . . . ,v m ) in the class II 9 . 

Proof: We first prove existence. Let A G R and for i G J , let Fi be the function defined by: 

Fi{t, x,y\...,y m ,z) = e xt fi(t, x, e~ xt yi , .., e~ xt y m , e~ xt Zi ) - X Vi . (5.14) 

Since fi is uniformly Lipschitz w.r.t. yi then Fi is so and for A small enough (A < 0) the function Fi 
is non-decreasing in all variables (y%, ...,y m ). Next thanks to Theorems Q] and [2l there exist processes 
(Y l , Z l , K l )i£j solution of the system of reflected BSDEs with interconnected obstacles associated 
with {(Fi(t,x,y l , ...,y m , z)) ieJ , (e XT hi(X^ x )) ieJ , (e Xt gij)ij e j). Additionally thanks to Proposition[37El 
there exist deterministic Isc functions v % , i G J , such that: 

Vi £j,Vt< T,Vs G [t,T\,Yj = v^X**). (5.15) 

Let us now analyze the decreasing scheme. First let us consider (Y, Z) the solution of the following 
standard BSDE: 

Y G S 2 , Z G V 24 

Y s = max i= i im e XT hi{X^ x ) + f g [max l= i )m F;](r, X l r ' x , Y r , . . . , Y r , Z r )dr - jj Z r dB r s <T. 

(5.16) 

Next for any i G JT, let us set Y h ° = Y and for n > 1 let us define (Y hn , Z hn , K hn ) by: 

' y«,n K i,n £ ^ z i,n £ n 2,d and R i, n non -decreasing; 

Yt n = e^hiiX^) + jf Fi(r, X? , (Y^ 1 )^, Z^ n )dr + K l > - K l s > n - jf Zp n dB r ; 

V:-" > max jeJ {Yr - < x *y, i {».Xl')}: " ' ' 

k fi(YP n - max ieJ {YJ-" - e Xs gi 3 (s,X^ x )})dKT = 0. 
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First note that the existence of Y is obvious by Pardoux-Peng's result [23] and then we easily deduce 
by induction: 

(i) for any n > 1, there exists a unique m-uplet of processes (Y l,n , Z t,n , K l ' n ) , i G J , ; 

(ii) for any n > 1, there exist deterministic continuous functions v l,n , i G J , such that: 

Vi G J, Vt < T,Vs G [t,T],y;- n = u*> n (s,^ s ); (5.18) 

(iii) for any i G J and n > 0, y*> n+1 < Y*>™ and u*> n > t;'' n+1 ; 

(iv) for any i £ J, the sequence (Y 4,n ) n converges to Y* in % 2,1 , where such a process has been 
introduced in equation (|5.15p . 

Actually for (i), we just need to use Theorem [1] and [2] since Y l, ° exists. To show (ii), we use induction: 
obviously the property is valid for n = (see e.g. [13], Thm.4.1). Next if the property holds for some 
n then it holds also for n + 1 in using Corollary [TJ As for (iii) , the property is true for n = since 
((y* , Z % , K % ) = (Y, Z, Q))i£j is the unique solution of the system associated with 

(fi := [maxi=i >m Fi\(r,Xr' x ,y l , ...,y m ,z)) ie j,(hi := max i= i jm e XT hi(X^ x )) i( zj, (e xt gij) i)je j) and then 
it is just enough to use the comparison result of Remark [1] Next if the property is valid for some n then 
it is also valid for n + 1, in using once more comparison since i*i(t, x, (y z )i£j, z), i G J , is non-decreasing 
in (y l )i & j. Finally (iv) holds true because the mapping $ defined in (|3.16p is a contraction in T-L 2 ' m 
and we obviously have: 

Vn > 0, (Y^+Vx,™ = <S>((Y*> n ) l=hm ). 

Let us now show that the deterministic functions (v l )i^j of (]5.15p are continuous. Recall that 
thanks to Proposition 13.11 each v % is the limit of some increasing sequence and hence it is Isc: therefore, 
it is enough to show that they are upper semicontinuous. But from the estimate given in (I3.17p . we 
know that: 

V(t,x)G [0,T] xM fc , \v%x)-v i > n {t,x)\=E[\Yt-Y*> n \] < C\\(Y% eJ - (Y^- 1 )^^. 

As limn^oc ||(y j )j ei 7 — (y i,ra_1 )j e: 7|| W 2,m = then for any i G J, the sequence {v l ' n ) n defined in (|5.17p 
converges pointwisely and decreasingly (from (|5.17p -(iii)) to v l . As a decreasing limit of continuous 
functions, v % is use and then continuous. It follows that the solution (y l ,z l ,k l )i ( zj of the system of 
reflected BSDEs with oblique associated with 

((fi(s,X s ' x ,y 1 , . . . ,y m ,z)) i( zj, (hi(X^ x )) ie j, (gij(s,X s ,x )) i)je j) has the following representation: 

Vi G J,Vt < T,Vs G [t,T\,Vx G R k , y\ = e" As ^(s, X l f). 

As v l , i G J , is continuous and of polynomial growth then using the result by El-Karoui et al.([12j, 
Thm. 8.5) related to connection between solutions of reflected BSDEs and viscosity solutions of PDE 
with obstacles we deduce that (e~ Ai f*(t, x))i & j is a viscosity solution for system f)2.3f) . ■ 
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As previously mentioned, we now consider the case when the functions —fi, i G S, verify (H2)-(iv). 
Then we have: 



Theorem 7 Assume that assumptions (H3), (H4) are fulfilled and that the functions (—fi)i^j verify 
(H2). Then the system of variational inequalities with inter- connected obstacles \2. 3\) has a continuous 
solution (y , . . . , v m ) in the class IT 9 . 

Proof: We first prove existence of a candidate to be a viscosity solution of the system. As previously, 
we will relate it to the unique solution of the multidimensional reflected BSDE (|3,3p . 

Step 1: Construction 

For i G J , let Fi be defined as in (|5.14p . Choosing A large enough, we obtain that each F l is 
decreasing with respect to all the variables yj, j = 1, • • • ,m. Next let us consider the following iterative 
Picard scheme: for any i G J", Y l, ° = and for n > 1, we define (Y %,n )i^j by: 

Therefore we know from Theorem [2] that the sequence (Y l,n )i^j converges in H 2 '" 1 to the unique 
solution (Y l )i e j of the system of reflected BSDEs with oblique reflection associated with 

Next using Proposition 13.11 and an induction argument, we deduce the existence of continuous 
functions with polynomial growth such that: 

Vn> 0,Vi G J,V(t,x) G [0, T] x R fc ,Vs G [t,T], Yp 71 = t/> n (s, X l / X ). 

But from Theorem [3}- (a) we have: for any i,n,m and s <T, 

E [\Y*' n - Yp m \ 2 ] < C\\(Y^) ie j - (Y^-^jW^ and 



E[|Y? - Y*' m \ 2 } < C\\(Y% eJ - (Y^) ie jf H%m , 



(5.19) 



where, to obtain the last inequality, we rely on the characterization of the solution (Y l )i^j constructed 
in [7J, i.e., (Y^j = $((Y*) i&7 ). Taking s = t we obtain: V (t,x) G [0,T] x R k , 

E[]Y i,n _ yi, m]] = lv i,n^ x) _ v i,m^ x) \ < C \\ (Y^) ieJ - (Y^^jW^. 

As the sequence ((Y l,n )i & j) n >Q is convergent in 1-L 2 ' m then it is of Cauchy type which implies that 
(V' n ) n >o is so and then converges pointwisely to a deterministic function v l , for any i G J . Thus going 
back to (|5.19p we deduce that: 

Vie J, VsG [t,T],¥-a.s., Y* = v i {s 1 X t *). (5.20) 

Let us now show that v l , i G J , belongs to IT 9 . Actually since $ is a contraction in ("H 2 ' m , ||.||/3 ) and 
by some induction procedure on n we get: 

/-in 

Vn,g > 0, \\{Y^) ieJ - (Y^WIk < —^-\\( Y ^) l&J \y o ; 
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with C<j> (such that < C$ < 1) which is the contraction constant of the mapping $ (constant which is 
independent of (t,x)). As the norms ||.|| and \\-\\/3 are equivalent then there exists a constant C\ such 
that: 

Vn,q > 0, \\(Y^) ieJ - (Y^ejW < C^IK^'Wll- 

Taking now the limit as q goes to +oo and in view of (15.190 and (|5.20[) . if we then take s = t we deduce 
that: 

V(t,x) G [0,T]xl fc , ^(t.a;) x)|<C 2 ||(y < ' 1 WII- 

Finally one can check easily that )jg^7"||(t,a;) is of polynomial growth (since E[sup s <ji 

belongs to II 9 for any 7 > 0, see (|3.2p ) and since v l,n is so, then we deduce that v % is also of polynomial 
growth for any i G J. 

Step 2: Continuity of v l , i G J ■ 

We again rely on the convergence result of any sequence (Y l,n ) n constructed via the Picard iterative 
scheme. So let us initialize the scheme as follows: 

Vi G {I,--. ,m},Vs <T, = C (l + \X l f\ p ) , (5.21) 

with the constant C and the integer p given by the fact that v % is in II 9 . 

Next let {Y l,t,x )i={i,— ,m} be the unique solution of the multidimensional RBSDE associated with 
((F t ( S ,xt' x ,y\...,y m ,z)) ie j,(e XT h i (X^)) teJ ,(e xt g^ then we have Y^ x = «<( s , Xt' x ) 

for any s G [t,T] and i ^ J . By definition of y*' in 15.211 and the hypothesis on v\ we obtain 

P-o.s.,Vi G J,Vs G [t, T], y?'* ,:c < Y*'°. 

Let us now prove by induction that for any n > we have: 

V n G N, G J, Vs G [t, T], yp n+1 < Y l s < Yp 2n . (5.22) 

For n = 0, the inequality of the right-hand side is already true. Let consider the left-hand side. Recalling 
that by construction we have (Y 1 ' 1 )^ = <&((Y l '°)j gi 7) and (Y l )j 6 j- = ^((Y l )i & j), we get the following 
comparison result for the two drivers of the RBSDE satisfied by Y l and Y*' 1 

Fi{s, X^,Y S \ Y™, z*) > Fi(s, xI*,yw, • • • , y s m '°, z% 

since Fi is non-increasing w.r.t. all variables yj, j G J . Next using the comparison result (see Remark 
HJ and uniqueness of the solution of the system associated with 

Jtf*,!/ 1 , ...,2/ m ,z)) ieJ , (e^hiiX^))^, (e^g^X^))^) we deduce that Y*' 1 < Y 4 for any 
i G JT. Thus the property (I5.22D is valid for n = 0. Now if it is satisfied for some n and repeating the 
same argumentation, it also holds for n + 1, whence the claim. 
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Next and relying once more on the result obtained El-Karoui et al.(|12j. Thm. 8.5) let v l ' n , i G J 
and n > 0, be the deterministic continuous functions of II 9 such that: 

V(i,x) G [0,T] x M fc ,Vs G [t,r],^' n = tJ i » n (s ) Jf*» x ). 

First and as previously for any i G J , the sequence (v l ' n ) n >o is of Cauchy type and converges pointwisely 
to v l . Next the inequalities (|5.22p imply that: Vn > and i G S, 

~i,2n+l < < 5 *,2» 

which implies 

u* = lim / ^' 2n+1 = lim \ ^' 2n . 

n n 

Therefore for any i G JT, v l is both txsc and Isc and thus it is continuous. Next, since (Y l , Z l , K l )i £ j 
is the unique solution of the system of reflected BSDEs associated with the following triplet of datas 
((Fiis, Xt>*,y\ ...,y m , z)) ieJ , (e AT ^(^) W, (e M 9ij (s, X*r))i d ej) then 

((e~ xt Yt,e~ xt Zl,e~ xt dKl)t<T)i<=j is the solution of the sytem of reflected BSDEs associated with 
((fi)iEj,(hi)i£j,(gij)i t j£j) then using once more the result by El-Karoui et al. Thm.8.5) to 

deduce that [e t v i )i^j is a continuous with polynomial growth solution of the system of variational 
inequalities with inter-connected obstacles (|2.3p . The proof is now complete. ■ 

Next we deal with the issue of uniqueness of the solution of (|2.3I) in the general case. 

Theorem 8 Under Assumptions (H1),(H2),(H3) and (H4), the solution of the system of variational 
inequalities with inter- connected obstacles \2. 3\) is unique in the class (vi, v m ) of continuous func- 
tions which belong to II 9 . 

Proof . As usual it is enough to show that if (m, ...,u m ) (resp. (vi, ...,v m )) is a continuous subsolution 
(resp. supersolution) of (|2.3p such that i G J, belong to li 9 then m < v l , for all i G J '. Thus 

classically we have uniqueness of f)2.3j) . 

Step 1: We first assume the existence of a constant A < —m. max{ Cj,j = 1, m} (Cj is the Lipschitz 
constant of fj involved in [H2]-(ii)) such that for any i G J ', fi verifies: 
V t,x,yi, ...,yi-i,y i+1 , ...,y m ,y,y, \iy>y then 

fi(t,x,yx, ...,yi-i,y,y i+ x, ...,y m ) - fi(t,x,yi, ...,yi-i,y,yi+i,...,y m ) < X(y - y). (5-23) 

Let 7 > and C be such that for any i G J we have: 

\ui(t,x)\ + ar) | < C(l + |x| 7 ), V(t,x) G [0,T] x R k . 

Next as in Lemma B~T1 for 9 > and v large enough (v l (t,x) + #e -z/ *|x| 27+2 )i = i jm is also a super- 
solution for ()3.3|) . Therefore it is enough to show that for any i G J7, we have: 

V(t,z) G [0,T] x R k , Ul {t,x) < v i (t,x) + 6e~ ut \x\ 2 "' +2 , 
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and taking the limit as — > we obtain the desired result. So let us set w' l)6,u {t,x) = v l (t,x) + 
6e^ ut \x\ 2l+2 ^ (t,x) G [0, T] x M fc which is denoted by w l for simplicity. Assume now that there exists a 
point (i,x) G [0,T] x M. k such that maxj £ j-(u*(t, x) — x)) > 0. Using the growth condition on Uj 
and w ! , there exists R > such that: 

Vi G J", V(t,x) G [0,T] x M fc s.t. |x| > R, Ui(t,x) -w\t,x) < 0. 

Taking into account the values of the subsolution and the supersolution at T, it implies that 

< max (t :E)e[0 T]x]Rfc max iG j ( Ui (t,x) - w l (t,x)) = 

max (i,x)e[o,T[xB(o,fi) max ieJ (ui(t,x) -w % (t,x)) = maxi^j (v,i(t* , x*) -w l (t*,x*)) 

where (t*,x*) G [0, T[xB(0, R). Now let J be as in (JO) and let j G J be such that 

Uj(t*,x*)> max (uk(t* , x*) — gjk(t* , x*)). (5.25) 

Next let $n(*j y) and i^ n , n > 0, be the same functions as in (|4,6p and let (t n , x n ,y n ) G [0, T]xB'(0, R) 2 
be the triple such that: 

<5>l(t n ,x n ,y n ) = max & n (t,x,y). 

(t,x,y)e[0,T]xB'(0,R) 2 

As in the proof of Theorem 14. 1[ one can show that: 

(t n , x n ,y n ) — > (t*,x*,x*) and n\x n — y n \ 2l+2 — > as n — > oo. (5.26) 

Now as Uj and are continuous functions and taking into account (|5.25p we deduce that for n large 
enough: 

Uj(t n ,x n )> max (u k (t n ,x n ) - gjk(t n ,x n )). (5.27) 

We next apply Crandall-Ishii-Lions's Lemma (see e.g. [8] or [13], pp.216) with <E>^,, Uj, and V9 n (recall 
that (|5.27p is satisfied) at (t n ,x n ,y n ), there exist (p",g",M™) G J 2 ' + (uj)(t n , x n ) and (p",<?™,M™) G 
J 2 ~{w^){t n ,y n ) such that: 

= d t (p n {tn,x n ,y n ) = 2(t n -t*), (resp. g") = d x cp n (t n , x n , y n )(resp.-d y (p n (t n , x n , y n )) and 

M »" ° )<A, + f4; (5.28) 

o -a™ y 2ti 

where = D 2 x y ^ip n {t n ,x n ,y n ). Taking now into account that {uiji^j (resp. (w l )i & j) is a subsolution 
(resp. supersolution) of (|2.3p and once more (|5.27p we deduce: 

-p£ - 6(t„,x n ) T .g™ - irr[(crCT T )(t ri ,x n )M^] - fj(t n ,x n , (ui(t n ,x n )) ie j,a(t n ,x n ) T .q£) < 



and 



-K-K*n,yn) T <-7:Tr[(aa T )(i n) :^ > 0. 



33 



Making the difference between those two inequalities yields: 

-{Pu-Pl) - (b(t n ,x n ) T .qZ - b(t n ,y n ) T .qZ) - ±Tr[{aa J (t n ,x n )M% - aa J (t n ,y n )M%}} 

-{fj(t n ,x n , (u i (t n ,x n )) i( zj,a(t n ,x n )' T .<?") - fj(t n ,y n , (w % (t n ,y n )) i( zj,a(t n ,y n ) T .<?")} < 

(5.29) 

and then 

ni X n ) .q u 

) - fj(t n ,x n , (w l (t n ,y n )) i£j ,a(t 

n i X n ) ■<?")) < A n 

where as in (|4,12p and (|4.13p we have lim n _ 5 . 00 A n < 0. Next linearizing fj, which is Lipschitz w.r.t. 
(y l )iej, and using (|5.23p we obtain: 

-X(uj(t n ,x n ) -w j (t n ,y n )) - ^ ®n k ((, u k{tm%n) -w k (t n ,y n )) < A n 

with <dil k which stands for the increment rate of fj w.r.t. y\. and which, thanks to the monotonicity 
assumption (see [H2] (iv)), is nonnegative and bounded by Cj the Lipschitz constant of fj. Thus 

-X(uj(t n ,x n ) - wi(t n ,y n )) < '(( u k(t n , x n ) - w k (t n ,y n )) + + A n 

< C fYlk&J-i((. U k( t n,X n ) -w k {t n ,y n )) + + A n . 

Taking now the limit as n — > oo and since j G J to obtain: 

- \( Uj (t* ,x*)- uP (t*,x*)) < C) Zkej-i K (f ,x*)-w k (t*,x*))+ 

< (m- l)Cj( Uj (t*,x*) - W i{t*,x*)) 

which contradictory since Uj(t*,x*) — w J (t*,x*) > and —A > mCj. Thus for any i G J , Ui < uf" . 

Step 2: The general case. 

For arbitrary A G R, if (uj)j^j (resp. (v^)j^j) is a subsolution (resp. supersolution) of (|2,3p then 
Uj(t,x) = e xt Uj(t,x) and v^(t,x) = e xt v^(t,x) is a subsolution (resp. supersolution) of the following 
system of variational inequalities with oblique reflection (see the proof of Proposition 14.14 Step 2): 

Viej, 

min < Vi(t, x) — max (—e xt gij(t, x) + Vj(t, x)) , 

-d t Vi(t, x) - Cvi(t, x) + Xvi(t, x) - e xt fi(t, x, {e~ xt v % {t, x)) ieJ , e- xt a T (t, x).D x v% x))} = 0; 
Vi(T,x) = e XT hi{x). 

(5.30) 

But in choosing A small enough the functions Fi(t,x, (ui)i^j,z) = —Xu + e xt fi(t,x, (e~ xt Ui)i^j,e~ xt z), 
i G J, satisfy condition (15.23H and then thanks to the result stated in Step 1, we have Ui < di, i G J . 
Thus Ui < Vi for any i G J which is the desired result. ■ 

As a by-product we have: 

Corollary 2 Assume that (HI), (H3), (H4) and (H2) are fulfilled, then there is a unique solution of 
the system \2. 3\) in the class of continuous functions with polynomial growth. ■ 
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